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Abstract 

In high-energy heavy-ion collisions, the correlations between the emitted 
particles can be used as a probe to gain insight into the charge creation 
mechanisms. In this article, we report the first results of such studies using 
the electric charge balance function in the relative pseudorapidity (A77) and 
azimuthal angle (Aip) in Pb-Pb collisions at v /J/\T/v — 2.76 TeV with the 
ALICE detector at the Large Hadron Collider. The width of the balance 
function decreases with growing centrality (i.e. for more central collisions) in 
both projections. This centrality dependence is not reproduced by HIJING, 
while AMPT, a model which incorporates strings and parton rescattering, 
exhibits qualitative agreement with the measured correlations in Atp but fails 
to describe the correlations in Arj. A thermal blast wave model incorporat- 
ing local charge conservation and tuned to describe the px spectra and v 2 
measurements reported by ALICE, is used to fit the centrality dependence 
of the width of the balance function and to extract the average separation 
of balancing charges at freeze-out. The comparison of our results with mea- 
surements at lower energies reveals an ordering with ^Snn'- the balance 
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functions become narrower with increasing energy for all centralities. This 
is consistent with the effect of larger radial flow at the LHC energies but 
also with the late stage creation scenario of balancing charges. However, the 
relative decrease of the balance function widths in At/ and Atp with central- 
ity from the highest SPS to the LHC energy exhibits only small differences. 
This observation cannot be interpreted solely within the framework where 
the majority of the charge is produced at a later stage in the evolution of the 
heavy-ion collision. 

Keywords: Balance function, charge correlations, ALICE LHC 



1. Introduction 

According to Quantum ChromoDynamics (QCD), the theory that de- 
scribes the strong interaction, at sufficiently high energy densities and tem- 
peratures, a new phase of matter exists in which the constituents, the quarks 
and the gluons, are deconfined jlj. This new state of matter is called the 
Quark Gluon Plasma (QGP). Its creation in the laboratory, the correspond- 
ing verification of its existence and the subsequent study of its properties are 
the main goals of the ultrarelativistic heavy-ion collision programs. Con- 
vincing experimental evidences for the existence of a deconfined phase have 
been published already at RHIC energies Q. Recently, the first experimental 
results from the heavy-ion program of the LHC experiments provided addi- 
tional indication (3), 0| for the existence of this state of matter at this new 
energy regime. 

Among the different observables, such as the anisotropic flow Q or the 
energy loss of high transverse momentum particles (H], the charge balance 
functions are suggested to be sensitive probes of the properties of the system, 
providing valuable insight into the charge creation mechanism and can be 
used to address fundamental questions concerning hadronization in heavy- 
ion collisions 

The system that is produced in a heavy-ion collision undergoes an ex- 
pansion, during which it exhibits collective behavior and can be described 
in terms of hydrodynamics 0|. A pair of particles of opposite charge that is 
created during this stage is subject to the collective motion of the system, 
which transforms the correlations in coordinate space into correlations in mo- 
mentum space. The subsequent rescattering phase after the hadronization 
will also affect the final measured degree of correlation. The balance function 
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being a sensitive probe of the balancing charge distribution in momentum 
space, quantifies these effects. The final degree of correlation is reflected in 
the balance function distribution and consequently in its width. It was sug- 
gested in ji| that narrow distributions correspond to a system that consists 
of particles that are created close to the end of the evolution. It was also 
suggested that a larger width may signal the creation of balancing charges 
at the first stages of the system's evolution 

The balance function describes the conditional probability that a particle 
in a bin Pi in momentum space will be accompanied (balanced) by a particle 
of opposite charge with momentum Pi. The general definition is given in 
Eq.H] 

B ab (P 2 \Pi) = \(c ab {P 2 \Pi) + C 6a (P 2 |Pi) 



-C bb {P 2 \P x )-C aa {P2\Pi)), (1) 

where C ab (p2\Pi) = N ab (P2\Pi) / N b (Pi) is the distribution of pairs of parti- 
cles, of type a and b, with momenta P2 and Pi, respectively, normalized to 
the number of particles b. Particles a and b could come from different particle 
species (e.g. 7r + -7r - , K + -K~, p-p). In this article, a refers to all positive and 
b to all negative particles. This analysis is performed for both particles in 
the pseudorapiity intervals \r)\ < 0.8. We assume that the balance function is 
invariant over pseudorapidity in this region, and report the results in terms 
of the relative pseudorapidity A77 = r) b — r) a and the relative azimuthal angle 
Acp = <p b — ip a , by averaging the balance function over the position of one of 
the particles (similar equation is used for B(A<p)): 



B + 4A V ) = 1(C+-(^V) + C-+(At7) - C__(At7) - C ++ (A V )). (2) 

Each term of Eq. [21 is corrected for detector and tracking inefficiencies as 
well as for acceptance effects, similar to what is proposed in Q, and can be 
written as C ab = (N ab /N b )/ f ab . The factors f ab (where in the case of charged 

particles, a and b correspond to the charge i.e. /_+, / ++ and / ) 

represent the probability that given a particle a is reconstructed, a second 
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particle emitted at a relative pseudorapidity or azimuthal angle (Arj or A(p, 
respectively), would also be detected. These terms are defined as the product 
of the single particle tracking efficiency e(r],(p,p'r) and the acceptance term 
a(Arj, Atp). The way they are extracted in this analysis is described in one 
of the following sections. 

For a neutral system, every charge has an opposite balancing partner and 
the balance function would integrate to unity. However, this normalization 
does not hold if not all charged particles are included in the calculation due 
to specific momentum range or particle type selection. 

The width of the balance function distribution can be used to quantify 
how tightly the balancing charges are correlated. It can be characterized 
by the average (Arj) or (A(p) in case of studies in pseudorapidity or the 
azimuthal angle, respectively. The mathematical expression for the case of 
correlations in pseudorapidity is given in Eq. [3] (similar for (Atp)). 



(3) 



i=i 



where (A^) is the balance function value for each bin Ar/j, with the sum 

running over all bins k. 

Experimentally, the balance function for non-identified particles was stud- 
ied by the STAR collaboration in Au-Au collisions at ^s^n = 130 GeV (sj( , 
followed by the NA49 experiment in Pb-Pb collisions at the highest SPS 
energy [sj. Both experiments reported the narrowing of the balance func- 
tion in At] in more central compared to peripheral collisions. The results 
were qualitatively in agreement with theoretical expectations for a system 
with a long-lived QGP phase and exhibiting delayed hadronization. These 
results triggere d an intense theoretical investigation of their interpretation 



10, 11, 12, 13, 14, 15, 16 



In [10J, it was suggested that the balance func- 
tion could be distorted by the excess of positive charges due to the protons 
of the incoming beams (unbalanced charges). This effect is expected to be 
reduced at higher collision energy, leaving a system at mid-rapidity that is 
net-baryon free. Also in [Ioj|, it was proposed to perform balance function 
studies in terms of the relative invariant momentum of the particle pair, to 



eliminate the sensitivity to collective flow. In [ill ], it was shown that purely 



hadronic models predict a modest broadening of the balance function for 
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central heavy-ion collisions, contrary to the experimentally measured nar- 
rowing. It was also shown that thermal models were in agreement with the 
(at that time) published data, concluding that charge conservation is local at 
freeze-out, consistent with the delayed charged-creation scenario fll| . Sim- 
ilar agreement with the STAR data was reported in [l2j|, where a thermal 
model that included resonances was used. In fl3jj , the author showed that 
the balance function, when measured in terms of the relative azimuthal an- 
gle of the pair, is a sensitive probe of the system's collective motion and 



in particular of its radial flow. In |l4j . it was suggested that radial flow is 



also the driving force of the narrowing of the balance function in pseudo- 
rapidity, with its width being inversely proportional to the transverse mass, 
mx = a/ m 2 + p\. In parallel in 15|, LL6| , the authors attributed the narrowing 
of the balance function for more central collisions to short range correlations 
in the QGP at freeze-out. 

Recently, the STAR collaboration extended their balance function studies 
in Au-Au collisions at ^/snn = 200 GeV fl7| , confirming the strong centrality 
dependence of the width in Arj but also revealing a similar dependence in A(p, 
the latter being mainly attributed to radial flow. Finally, in [l8| the authors 
fitted the experimentally measured balance function at the top RHIC energies 
with a blast -wave parameterization and argued that in Atp the results could 
be explained by larger radial flow in more central collisions. However the 
results in At] could only be reproduced when considering the separation of 
charges at freeze-out implemented in the model. They also stressed the 
importance of performing a multi-dimensional analysis. In particular, they 
presented how the balance function measured with respect to the orientation 
of the reaction plane (i.e. the plane of symmetry of a collision defined by the 
impact parameter vector and the beam direction) could probe potentially 
one of the largest sources of b ackg round in studies related to parity violating 
effects in heavy-ion collisions (l9j. 

In this article we report the first results of the balance function measure- 
ments in Pb-Pb collisions at y / J/v r /v — 2.76 TeV with the ALICE detector 



20l [21] . The article is organized as follows: Section [5] briefly describes the 



experimental setup, while details about the data analysis are presented in 
Section [3j In Section H] we discuss the main results followed by a detailed 
comparison with different models in Section [5) In the same section we present 
the energy dependence of the balance function. We conclude with the sum- 
mary and a short outlook. 
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2. Experimental setup 



ALICE [21| is the dedicated heavy-ion detector at the LHC, designed 
to cope with the high charged-particle densities measured in central Pb- 
Pb collisions j22|. The experiment consists of a large number of detector 
subsystems inside a solenoidal magnet (0.5 T). The central tracking systems 
of ALICE provide full azimuthal coverage within a pseudorapidity window 
1 7/| < 0.9. They are also optimized to provide good momentum resolution 
(~ 1% at Pt < 1 GeV/c) and particle identification (PID) over a broad 
momentum range, the latter being important for the future, particle type 
dependent balance function studies. 

For this analysis, the charged particles were reconstructed using the Time 
Projection Chamber (TPC) [23], which is the main tracking detector of the 
central barrel. In addition, a complementary analysis relying on the combined 
tracking of the TPC and the Inner Tracking System (ITS) was performed. 
The ITS consists of six layers of silicon detectors employing three different 
technologies. The two innermost layers are Silicon Pixel Detectors (SPD), 
followed by two layers of Silicon Drift Detectors (SDD). Finally the two 
outermost layers are double-sided Silicon Strip Detectors (SSD). 

The position of the primary interaction was determined by the TPC and 
by the SPD, depending on the tracking mode used. A set of forward detec- 
tors, namely the VZERO scintillator arrays, were used in the trigger logic 
and also for the centrality determination [22|. The VZERO detector consists 
of two arrays of scintillator counters, the VZERO-A and the VZERO-C, po- 
sitioned on each side of the interaction point. They cover the pseudorapidity 
ranges of 2.8 < rj < 5.1 and -3.7 < 77 < -1.7 for VZERO-A and VZERO-C, 
respectively. 



For more details on the ALICE experimental setup, see 
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3. Data analysis 

Approximately 15 x 10 6 Pb-Pb events, recorded during the first LHC 
heavy-ion run in 2010 at ^Snn = 2.76 TeV, were analyzed. A minimum bias 
trigger was used, requiring two pixel chips hit in the SPD in coincidence with 
a signal in the VZERO-A and VZERO-C detectors. Measurements were also 
made with the requirement changed to a coincidence between signals from 
the two sides of the VZERO detectors. An offline event selection was also 
applied in order to reduce the contamination from background events, such as 
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electromagnetic and beam-gas interactions. All events were required to have 
a reconstructed vertex position along the beam axis (V z ) with \V Z \ < 10 cm 
from the nominal interaction point. 




Figure 1: (color online). The correction factor /_| (A77, A(p) for the 5% most central Pb- 

Pb collisions, extracted from the single particle tracking efficiencies e(r),ip,pr) and the 
acceptance terms a(Arj, Atp) (see text for details). 



The data were sorted according to centrality classes, reflecting the ge- 
ometry of the collision (i.e. impact parameter), which span — 80% of the 
inelastic cross section. The — 5% bin corresponds to the most central 
(i.e. small impact parameter) and the 70 — 80% class to the most periph- 
eral (i.e. large impact parameter) collisions. The centrality of the collision 
was estimated using the charged particle multiplicity distribution and the 
distribution of signals from the VZERO scintillator detectors. Fitting these 



distributions with a Glauber model [241 ]. the centrality classes are mapped to 
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Figure 2: (color online). Balance function as a function of Ar] for different centrality 
classes: 0-5% (a), 30-40% (b) and 70-80%. Mixed events results, not corrected for the 
detector effects, are shown by open squares. See text for details. 



the corresponding mean number of participating nucleons (N part ) [25] . Dif- 
ferent centrality estimators (i.e. TPC tracks, SPD clusters) were used to 
investigate the systematic uncertainties. 

To select charged particles with high efficiency and to minimize the con- 
tribution from background tracks (i.e. secondary particles originating either 
from weak decays or from the interaction of particles with the material), all 
selected tracks were required to have at least 70 reconstructed space points 
out of the maximum of 159 possible in the TPC. The (x 2 ) per degree of 
freedom the momentum fit was required to be below 2. To further reduce 
the contamination from background tracks, a cut on the distance of clos- 
est approach between the tracks and the primary vertex (dca) was applied 
(dca xy /d xy ) 2 + (dca z /d z ) 2 < 1 with d xy = 2.4 cm and d z = 3.2 cm. In the 
parallel analysis, with the combined tracking of the TPC and the ITS, the 
values of d xy = 0.3 cm and d z = 0.3 cm were used, profiting from the better 
dca resolution that the ITS provides. Finally, we report the results for the 
region of |?7| < 0.8 and 0.3 < Pt < 1.5 GeV/c. The p^ range is chosen to 
ensure a high tracking efficiency (lower cut) and a minimum contribution 
from (mini-)jet correlations (upper cut). 
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Figure 3: (color online). Balance function as a function of Aip for different centrality 
classes: 0-5% (a), 30-40% (b) and 70-80%. Mixed events results, not corrected for the 
detector effects, are shown by open squares. See text for details. 



4. Results 

As discussed in the introduction, the correction factors /+_, / |_, /++, 

and / are needed to eliminate the dependence of the balance function on 

the detector acceptance and tracking inefficiencies. The tracking inefficiencies 
are extracted from a detailed Monte Carlo simulation of the ALICE detec- 
tor based on GEANT3 j2fjj|. The acceptance part of the correction factors, 
en (At/, Aip), is extracted from mixed events. The mixed events are generated 
by taking all two-particle non-same-event combinations for a collection of 
a few (« 5) events with similar values of the z position of the reconstructed 
vertex (|A14| < 5 cm). In addition, the events used for the event mixing be- 
longed to the same centrality class and had multiplicities that did not differ 
by more than 1-2%, depending on the centrality. Figure [1] presents the cor- 
rection factor for the distribution of pairs of particles with opposite charge as 
a function of the relative pseudorapidity and azimuthal angle differences for 
the 5% most central Pb-Pb collisions. The maximum value is observed for 
At/ = and is equal to the px^integrated single particle efficiency. The distri- 
bution decreases to ps near the edge of the acceptance i.e. \Ar]\ m 1.6. This 
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reduction reflects the decrease of the probability of detecting both balancing 
charges as the relative pseudorapidity difference increases. The correction 
factor is constant as a function of Atp. 

The measured balance function is averaged over positive and negative 
values of Arj (Atp) and reported only for positive values. The integrals of the 
balance function over the reported region are close to 0.5, reflecting the fact 
that most of the balancing charges are distributed in the measured region. 

Figure [2] presents the balance functions as a function of the relative pseu- 
dorapidity Arj for three different centrality classes: the 0-5% (most central), 
the 30-40% (mid-central) and the 70-80% (most peripheral) centrality bins. 
It is seen that the balance function, in full circles, gets narrower for more cen- 
tral collisions. Figure [2] presents also the balance functions for mixed events, 
not corrected for detector effects, represented by the open squares. These bal- 
ance functions, fluctuate around zero as expected for a totally uncorrelated 
sample where the charge is not conserved. 

Figure [3] presents the balance functions as a function of the relative az- 
imuthal angle for the same centrality classes as in Fig. |5J The balance func- 
tions calculated using mixed events and not corrected for the tracking effi- 
ciency exhibit a distinct modulation originating from the 18 sectors of the 
TPC. This modulation is more pronounced for more central collisions, since 
the charge dependent acceptance differences scale with multiplicity. The 
efficiency-corrected balance functions, represented by the full markers, indi- 
cate that these detector effects are successfully removed. Narrowing of the 
balance function in more central events has been also observed in this repre- 
sentation. A decrease of the balance function at small Aip (i.e. for Atp < 10°) 
can be observed for the mid-central and peripheral collisions. This can be 
attributed to short-range correlations between pairs of same and opposite 
charge, such as HBT and Coulomb effects fli| . 

In both Fig. [2] and Fig. [3] as well as in the next figures, the error bar 
of each point corresponds to the statistical uncertainty (typically the size of 
the marker). The systematic uncertainty is represented by the shaded band 
around each point. The origin and the value of the assigned systematic un- 
certainty on the width of the balance function, calculated for each centrality 
and for both At] and Aip, will be discussed in the next paragraph. 

The data sample was analyzed separately for two magnetic field config- 
urations. The two data samples had comparable statistics. The maximum 
value of the systematic uncertainty, defined as half of the difference between 
the balance functions in these two cases, is found to be less than 1.3% over 



21 



all centralities. In addition, we estimated the contribution to the system- 
atic uncertainty originating from the centrality selection, by determining 
the centrality not only with the VZERO detector but alternatively using 
the multiplicity of the TPC tracks or the number of clusters of the second 
SPD layer. This resulted in an additional maximum contribution to the esti- 
mated systematic uncertainty of 0.8% over all centralities. Furthermore, we 
investigated the influence of the ranges of the cuts in parameters such as the 
position of the primary vertex in the z coordinate, the dca, and the number 
of required TPC clusters. This was done by varying the relevant ranges, 
one at a time, and again assigning half of the difference between the lower 
and higher value of the width to the systematic uncertainty. The maximum 
contribution from these sources was estimated to be 1.3%, 1.1% and 1.3% 
for the three parameters, respectively. We also studied the influence of the 
different tracking modes used by repeating the analysis using tracks recon- 
structed by the combination of the TPC and the ITS (global tracking). The 
resulting maximum contribution to the systematic uncertainty of the width 
from this source is 1.1%, again over all centralities. Finally, the applied 
acceptance corrections result in large fluctuations of the balance function 
points for some centralities towards the edge of the acceptance (i.e. large 
values of Arj), which originates from the division of two small numbers. To 
account for this, we average over several bins at these high values of At] to 
extract the weighted average. This procedure results in an uncertainty that 
has a maximum value of 5% over all centralities. All these contributions are 
summarized in Table [TJ The final systematic uncertainty for each centrality 
bin was calculated by adding all the different sources in quadrature. The 
resulting values for the 0-5%, 30-40% and 70-80% centrality bins were es- 
timated to be 2.5%, 3.0% and 3.6%, respectively, in (Arj) (1.9%, 1.2% and 
2.4%, respectively, in (Acp)). 

5. Discussion 

5.1. Centrality dependence 

The width of the balance function (Eq. [3]) as a function of the centrality 
percentile is presented in Fig. HI Central (peripheral) collisions correspond 
to small (large) centrality percentile. The width is calculated in the entire 
interval where the balance function was measured (i.e. 0.0 < At] < 1.6 
and 0° < Aip < 180°). Both results in terms of correlations in the relative 
pseudorapidity ((Ar/)-upper panel, Fig.|4]-a) and the relative azimuthal angle 
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Table 1: The maximum value of the systematic uncertainties on the width of the balance 
function over all centralities for each of the sources studied. 



Systematic Uncertainty 


Category 


Source 


Value (max) 


Magnetic field 


(++)/(- -) 


1.3% 


Centrality estimator VZERO, TPC, SPD 


0.8% 




dca 


1.3% 


Cut variation 


N duster s\T PC) 


1.1% 




AV Z 


1.3% 


Tracking 


TPC, Global 


1.1% 


Binning 


Extrapolation to large Arj 


5.0% 



((A^-lower panel, Fig. H]-b) are shown. The experimental data points, 
represented by the full red circles, exhibit a strong centrality dependence: 
more central collisions correspond to narrower distributions (i.e. moving 
from right to left along the x-axis) for both At] and A<p. Our results are 
compared to different model predictions, such as HIJING |27| and different 
versions of a multi-phase transport model (AMPT) j28j. The error bars in 
the results from these models represent the statistical uncertainties. 

The points from the analysis of HIJING Pb-Pb events at ^s^n — 2.76 TeV, 
represented by the blue triangles, show little centrality dependence in both 
projections. The slightly narrower balance functions for central collisions 
might be related to the fact that HIJING is not just a simple superposition 
of single pp collisions; jet-like effects as well as increased resonance yields in 
central collisions could be reflected as additional correlations. The balance 
function widths generated by HIJING are much larger than those measured 
in the data, consistent with the fact that the model lacks collective flow. 

In addition, we compare our data points to the results from the analysis 
of events from three different versions of AMPT in Fig. HI The AMPT model 
consists of two different configurations: the default and the string melting. 
Both are based on HIJING to describe the initial conditions. The partonic 
evolution is described by the Zhang's parton cascade (ZPC) [29j|. In the de- 
fault AMPT model, partons are recombined with their parent strings when 
they stop interacting, and the resulting strings are converted to hadrons us- 
ing the Lund string fragmentation model. In the string melting configuration 
a quark coalescence model is used instead to combine partons into hadrons. 
The final part of the whole process, common between the two configura- 
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Figure 4: (color online). The centrality dependence of the width of the balance function 
(Arj) and (Aip), for the correlations studied in terms of the relative pseudorapidity and 
the relative azimuthal angle, respectively. The data points are compared to the predictions 
from HIJING [27:], and AMPT [H]. 
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Figure 5: (color online). The balance functions for the 5% most central Pb-Pb collisions 
measured by ALICE as a function of the relative pseudorapidity (a) and the relative 
azimuthal angle (b). The experimental points are compared to predictions from HIJING 
|27j . AMPT [2i| and from a thermal blast wave (3l|, 
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tions, consists of the hadronic rescattering which also includes the decay of 
resonances. 

The filled green squares represent the results of the analysis of the string 
melting AMPT events with parameters tuned j30[ to reproduce the measured 
elliptic flow (V2) values of non-identified particles at the LHC [jjj]. The width 
of the balance functions when studied in terms of the relative pseudorapidity 
exhibit little centrality dependence despite the fact that the produced system 
exhibits significant collective behavior 30]. However, the width of the balance 



function in Aip is in qualitative agreement with the centrality dependence of 
the experimental points. This is consistent with the expectation that the 
balance function when studied as a function of A09 can be used as a measure 



of radial flow of the system, as suggested in |13l . Il8| . We also studied the 
same AMPT configuration, i.e. the string melting, this time switching off 
the last part where the hadronic rescattering takes place, without altering 
the decay of resonances. The resulting points, indicated with the orange 
filled stars in Fig. HI demonstrate a similar qualitative behavior as in the 
previous case: no centrality dependence of (Arj) and a significant decrease 
of (Aip) for central collisions. On a quantitative level though, the widths in 
both projections are larger than the ones obtained in the case where hadronic 
rescattering is included. This can be explained by the fact that within this 
model, a significant part of radial flow of the system is built during this very 
last stage of the system's evolution. Therefore, the results are consistent 
with the picture of having the balancing charges more focused under the 
influence of this collective motion, which is reflected in a narrower balance 
function distribution. In addition, we analyzed AMPT events produced using 
the default configuration, which results in smaller v n flow coefficients but 
harder spectra than the string melting. The extracted widths of the balance 
functions are represented by the open brown squares and exhibit similar 
behavior as the results from the string melting configuration. In particular, 
the width in Ar) shows little centrality dependence while the values are in 
agreement with the ones calculated from the string melting. The width in 
A(p shows similar (within the statistical uncertainties) quantitative centrality 
dependence as the experimental data points. However, these latter results 
exhibit a systematically lower value of the width for this version of AMPT 
compared to the string melting configuration. This latter effect is consistent 
with the observation of having a system exhibiting larger radial flow with 
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the default versionH 

Finally, we fit the experimentally measured values with a thermal blast- 



wave model |3lL |32|. This model, assumes that the radial expansion velocity 
is proportional to the distance from the center of the system and takes into 
account the resonance production and decay It also incorporates the lo- 
cal charge conservation, by generating ensembles of particles with zero total 
charge. Each particle of an ensemble is emitted by a fluid element with a 
common collective velocity following the single-particle blast-wave parame- 
terization with the additional constraint of being emitted with a separation 
at kinetic freeze-out from the neighboring particle sampled from a Gaussian 
with a width denoted as a n and a v in the pseudorapidity space and the az- 
imuthal angle, respectively. The procedure that we followed started from 
tuning the input parameters of the model to match the average values 
extracted from the analysis of identified particle spectra 34| as well as the 
V2 values for non-identified particles reported by ALICE (J\Ve then adjust 
the widths of the parameters a v and a v to match the experimentally mea- 
sured widths of the balance function, (At/) and (Ay?) . The resulting values 
of a v and a v are listed in Table [2J We find that a rj starts from 0.28 ±0.05 for 
the most central Pb-Pb collisions reaching 0.52 ± 0.07 for the most periph- 
eral, while starts from 0.30 ± 0.10 evolving to 0.76 ± 0.01 for the 60-70% 
centrality bin. 

Figure presents the detailed comparison of the model results with the 
measured balance functions as a function of At/ (a) and Aip (b) for the 
5% most central Pb-Pb collisions. The data points are represented by the 
full markers and are compared with HIJING (dotted blue line), AMPT string 
melting (dashed green line) and the thermal blast-wave (full black line). The 
distributions for HIJING and AMPT are normalized to the same integral to 
facilitate the direct comparison of the shapes and the widths. It is seen 
that for correlations in the relative pseudorapidity, both HIJING and AMPT 
result in similarly wider distributions. As mentioned before, the blast-wave 
model is tuned to reproduce the experimental points, so it is not surprising 
that the relevant curve not only reproduces the same narrow distribution but 
describes fairly well also its shape. For the correlations in A(p the HIJING 



3 We recently confirmed that AMPT does not conserve the charge. The influence of 
this effect to our measurement can not be easily quantified. However we still consider 
interesting and worthwhile to point out that this model describes in a qualitative (and to 
some extent quantitative) way the centrality dependence of (Aip). 
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Table 2: 



The values of and cr v extracted by fitting the centrality dependence of both 



j i i I 

( A77) and (Aip) with the blast-wave parameterization of 3l|, |32j • 



Results from the fit with the blast-wave model 



Centrality 


a v 


<J 9 


0-5% 


0.28 ±0.05 


0.30 ±0.10 


5-10% 


0.32 ±0.05 


0.35 ±0.07 


10-20% 


0.31 ±0.05 


0.36 ±0.08 


20-30% 


0.36 ±0.03 


0.43 ±0.05 


30-40% 


0.43 ± 0.04 


0.52 ±0.05 


40-50% 


0.42 ± 0.04 


0.54 ±0.06 


50-60% 


0.44 ±0.07 


0.64 ±0.06 


60-70% 


0.52 ±0.07 


0.76 ±0.01 



curve clearly results in a wider balance function distribution. On the other 
hand, there is a very good agreement between the AMPT curve and the 
measured points, with the exception of the first bins (i.e. small relative 
azimuthal angles) where the magnitude of B + -(A<p) is significantly larger 
in real data. This suggests that there are additional correlations present in 
these small ranges of Aip in data than what the model predicts. 

5.2. Energy dependence 

Figure [6] presents the comparison of our results for the centrality depen- 
dence (i.e. as a function of the centrality percentile) of the width of the bal- 
ance function, (A77) (Fig. [6]-a) and (A(p) (Fig.[B]-b), with results from STAR 
(l7j in Au-Au collisions at ^/snn — 200 GeV (stars). The ALICE points 
have been corrected for acceptance and detector effects, using the correction 
factors f a b, discussed in the introduction. To make a proper comparison with 
the STAR measurement, where such a correction was not applied, we employ 
the procedure suggested in Q to the RHIC points. Based on the assumption 
of a boost-invariant system the balance function studied in a given pseudo- 
rapidity window 5 + _(A?7|?7 max ) can be related to the balance function for an 
infinite interval according to the formula of Eq. H] 



B + -(A*m») = B + 4A V \oo) ■ (l - -^-). (4) 
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This procedure results in similar corrections as to the case where the / a j are 
used, if the acceptance is flat in 77 (which is a reasonable assumption for the 
acceptance of STAR)0 

While the centrality dependence is similar for both measurements, the 
widths are seen to be significantly narrower at the LHC energies. This is 
consistent with the idea of having a system exhibiting larger radial flow at 
the LHC with respect to RHIC |3| while having a longer-lived QGP phase 
(33| with the consequence of a smaller separation between charge pairs when 
created at hadronization. However, it is seen that the relative decrease of the 
width between central and peripheral collisions seems to be similar between 
the two energies. This observation could challenge the interpretation of the 
narrowing of the width in A77 as primarily due to the late stage creation of 
balancing charges. 

To further quantify the previous observation, Fig. [7] presents the relative 
decrease of (At]) (a) and (A<p) (b) from peripheral to central collisions as a 
function of the mean number of participating nucleons, (iV part ), for the high- 
est j^l and RHIC |l7| energies, compared to the values reported in this 
article. In this figure, central (peripheral) collisions correspond to high (low) 
number of (N pa , Tt ). The choice of the representation as a function of (iVp art ) is 
mainly driven by the apparent better scaling compared to the centrality per- 
centile. It is seen that in terms of correlations in relative pseudorapidity the 
data points at the different energies are in fairly good agreement within the 
uncertainties, resulting though into an additional, marginal decrease for the 
0-5% most central collisions of m (9.5 ± 2.0 (stat) ± 2.5 (syst))% compared 
to the RHIC point. On the other hand, (Atp) j (Atp) peripheral exhibits a de- 
crease of w (14.0 ± 1.3 (stat) ± 1.9 (syst))% between the most central Au-Au 
collisions at ^/snn — 200 GeV and the results reported in this article. This 
could be attributed to the additional increase in radial flow between central 
and peripheral collisions at the LHC compared to RHIC energies. Another 
contribution might come from the bigger influence from jet-like structures at 
the LHC with respect to RHIC that results in particles being emitted prefer- 



4 We do not compare our results to the data from the NA49 experiment at SPS in 
this figure, for two reasons. Firstly, the balance function in that experiment was not 
measured at mid-rapidity. Secondly, the non-uniform acceptance in pseudorapidity makes 
the simplified correction of Eq. [4] invalid. 

5 We include the NA49 points in this representation since the ratio to the peripheral 
results should cancel out the acceptance effects to first order. 
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Figure 6: (color online). The centrality dependence of the balance function width (Arf) 
(a) and (Acp) (b). The ALICE points are compared to results from STAR 0. The STAR 
results have been corrected for the finite acceptance as suggested in Q. 
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Figure 7: (color online). The centrality dependence of the relative decrease of the width 
of the balance function in the relative pseudorapidity (a) and relative azimuthal angle (b). 
The ALICE points are compared to results for the highest SPS [j| and RHIC [ItJ energies. 
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entially in cones with small opening angles. Contrary to (Aip) / (Atp) peripheral, 
this strikingly marginal decrease of (A77) / {At]) peripheral between the three col- 
liding energy regimes that differ more than an order of magnitude, can not 
be easily understood solely within the framework of the late stage creation 
of charges. 

6. Summary 

This article reported the first measurements of the balance function for 
charged particles in Pb-Pb collisions at the LHC using the ALICE detector. 
The balance function was studied both, in relative pseudorapidity (At]) and 
azimuthal angle (Atp). The widths of the balance functions, (At]) and (Atp), 
are found to decrease when moving from peripheral to central collisions. The 
results are consistent with the picture of a system exhibiting larger radial 
flow in central collisions but also whose charges are created at a later stage 
of the collision. While HIJING is not able to reproduce the observed central- 
ity dependence of the width in either projection, AMPT tuned to describe 
the v 2 values reported by ALICE seems to agree qualitatively with the cen- 
trality dependence of (Atp) but fails to reproduce the dependence of (Ar]). 
A thermal blast-wave model incorporating the principle of local charge con- 
servation was fitted to the centrality dependence of (At]) and (Atp). The 
resulting values of the charge separation at freeze-out can be used to con- 
strain models describing the hadronization processes. The comparison of the 
results with those from lower energies showed that the centrality dependence 
of the width, in both the relative pseudorapidity and azimuthal angle, when 
scaled by the most peripheral widths, exhibits minor differences between 
RHIC and LHC. 

These studies will soon be complemented by and extended to the correla- 
tions of identified particles in an attempt to probe the chemical evolution of 
the produced system, to quantify the influence of radial flow to the narrow- 
ing of the balance function width in more central collisions and to further 
constrain the parameters of the models used to describe heavy-ion collisions. 
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